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We summarize recent results by the BaBab. and Belle Collaborations on the search for CP violation 
and new physics beyond the Standard Model in several rare B meson decays. The measurements 
exploit the final BaBar and Belle T(4S) dataset, which consist of approximately 470 and 772 million 
BB pairs, respectively. 

I. INTRODUCTION 

The B factory experiments, BABAR at SLAC and Belle at KEK, have been operated successfully in the first 
decade of the 21st century, collecting in total more than one billion pairs of charged {B+B-) and neutral (B°B ) B 
mesons. Such large samples have been obtained from the decays of boosted T(4S*) mesons produced in asymmetric 
e + e~ collisions at y/s = mx, e + e~ — > T(4S') — > BB. After establishing experimentally CP violation in the B 
meson system and confirming the Cabibbo-Kobayashi-Maskawa (CKM) paradigm [212] by measuring with good or 
excellent precision the sides and angles of the unitarity triangle in the "golden" modes, like B — > ccK^° for the 
angle (3 — arg[— V c dV* b /V t dV t * b ] (sin 2(3 — 0.679 ± 0.020 ;3j), a significant fraction of the B factories' physics program 
has been focused on rare B decays, which can provide additional constraints on the unitarity triangle parameters, or 
can probe the presence of new physics (NP) beyond the Standard Model (SM). 

II. EXPERIMENTAL TECHNIQUES 

Both BABAR and Belle are multi-purpose, (longitudinally) asymmetric detectors with nearly 4-7T coverage. In the 
e + e~ collisions the T(AS) are produced with a Lorentz boost "f/3 along the beam axis of 0.56 in BABAR and 0.42 
in Belle, resulting in an average longitudinal separation between the two B decay vertices of « 200-250 microns, 
from which the proper time difference between the two decays can be inferred. The two detectors consist of an 
inner silicon detector to reconstruct B and D decay vertices and low-p^ charged particles, a gaseous drift chamber 
for charged particle tracking, various devices optimized for excellent charged particle identification (PID) at low 
energy (w 1 GeV) through various techniques (dE/dx, Cherenkov light, time of flight, ..), homogeneous scintillating 
electromagnetic calorimeters for e/7 reconstruction and identification, and muon and Kl detectors installed in the 
iron yoke used to contain the flux return of the solenoid that provides the magnetic field for tracking. 

BB events produce typically a few number of charged (~ 11) and neutral particles and often it is possible to 
fully reconstruct a B decay, i.e. to measure directly each of its decay products. This allows the calculation of 
two kinematic quantities characterizing the B meson, the energy-substituted invariant mass (toes or Mf, c ), toes = 
\A-^t>cam) 2 — (p*b) 2 > an d the energy difference, AE = E* B — £^ cam - For correctly reconstructed B decays, these 
two quantities peak around the B meson mass (tob) and zero, respectively. Sometimes the B decay under study 
produces "invisible" particles that do not interact with the detector, for instance a neutrino or an hypothetical weakly 
interacting supersymmetric particle: in that case, one can still fully reconstruct the recoiling B (usually denoted as 
Stag) and use it to infer the four-momentum of the signal B and of the invisible particles. 

When a pair of neutral B mesons is produced, the flavors of the two mesons are entangled. If the first B meson 
decays at time t to a flavor specific final state, which can be identified through the presence of leptons, kaons or soft 
pions, the second B meson is inferred to have at the same time t the opposite flavor. These correlations are exploited 
in time-dependent CP violation measurements. 

The e + e~ collisions can also produce light quark-antiquark pairs, qq (q — u,d,s,c). The BB production cross 
section is around one fifth of the total hadronic cross section. The hadronic qq backgrounds are typically distinguished 
from BB events by means of event-shape variables, that discriminate between the "jet-like" topology of qq events 
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TABLE I: Results of the B — > nh measurements by the Belle Collaboration [13] . 



Decay 


BF (10 -6 ) 


BF significance (er) 


A C p 


Acp significance (a) 


B+ -> nK+ 


2.12 ±0.23 ±0.11 


13.2 


-0.38 ±0.11 ±0.01 


3.8 


B+ -> 777T + 


4.07 ±0.26 ±0.21 


22.4 


-0.19 ±0.06 ±0.01 


3.0 


B° ->■ r,^ 


1.27±g;|| ± 0.08 


5.4 







and the more "spherical" BB events, and also exploit the different angular distributions expected from angular 
momentum conservation and spin arguments. Because of the large correlations between the shape variables, those 
are usually combined in a single discriminating quantity using multivariate techniques (Fisher discriminants, neural 
networks), trained on clean samples of qq events collected at y/s « 40 — 50 MeV below the T(45) peak and with 
simulated signal samples. 

III. SEARCH FOR CP VIOLATION IN CHARMLESS B DECAYS 

Charmless hadronic B decays provide sensitive probes for CP violation in the Standard Model and beyond, since: 

• large CP-violating effects are predicted by the SM in some of these channels, and can therefore be measured 
with good precision; 

• the amount of CP violation in some charmless B decays is related, in the SM, to the amount of CP violation 
- measured with excellent precision - in B° — > (cc)K° decays; 

• some of these decays proceeds mainly (or purely) through loop-mediated diagrams and SM predictions could 
therefore be significantly affected by possible new physics contributions appearing in the loops of these diagrams. 

A. Search for CP violation in two-body B — > nh (h = tt + , K + , K°) decays 

Calculations of the SM contributions to the B — > r/K and B —> nir decays, combined with the measured values of 
the B^ — > rj(''K ± and B° — > n'K° branching fractions, lead to the following expectations [4Tll2j: 

1. the branching fraction of B° — > ijK is expected to be lower than that of B ± — > nK ± , because the tree diagram 
in the B° — » r/K° decay is color suppressed. 

2. in B ± — > r/K ± , a large direct CP asymmetry (Acp) is expected, as a consequence of the interference between 
the destructive combination of penguin amplitudes and the tree amplitude. 

3. in B ± — > rj7r , a large Aqp is also anticipated, due to the interference between the b — > d penguin and b — > u 
tree diagrams. 

Belle has published in [13] an updated measurement of B — > ijh (h — ir + , K + , K°) decays[36] using the final data 
sample of 772 x 10 6 BB pairs. The i] meson is reconstructed in both the 7r + 7r _ 7r° and 77 final states, while the K° 
meson is selected as a Kg decaying to two oppositely-charged pions. The full B decay tree is reconstructed, with 
an efficiency - including secondary branching fractions - of around 1.5% (4%) for the r/K and 5% (14%) for the 
r]h + final states, in the i] — > 7r + 7r~7r (7/ — > 77) channel. The signal yield is obtained through a maximum likelihood 
(ML) fit to Mfc c , AE and an event-shape variable R' . The samples containing the nK + and r)ii + candidates are 
fitted simultaneously, in order to constrain the background from the nK -f-> nn cross-feed. The projections of the 
fit on the Mf, c variable for B + — s- nK + and B~ — > i]K~ candidates, selected after further enhancing S/B with the 
requirements —0.1 < AE < 0.08 GeV and R' > 1.95, are shown in Figure [I] 
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FIG. 1: Mbc projections for B + — > nK + (left) and B~ — > nK~ (right) candidate events selected by Belle [13) . Points with 
error bars represent the data, the total fit functions are shown by black solid curves, signals are shown by red solid curves, 
dashed curves are the continuum contributions, dotted curves are the cross-feed backgrounds from misidentification and filled 
histograms are the contributions from charmless B backgrounds. The projections are for events that have —0.10 < AE < 0.08 
GeV and R' > 1.95. 



The main results are summarized in Table [T| Both B + — > t]K + and B + — » rjir + signals are observed with large 
statistical significances, exceeding 10a and 20a respectively, and both B decays exhibit large negative CP asymme- 
tries, around —40% and —20% respectively, at least 3a different from zero. The B — > nK° decay is observed for the 
first time, with a significance of more than 5a, and its branching fraction is smaller than that of B + — > nK + . All 
results are in agreement with SM expectations. 



B. Search for CP violation in three-body B — > KKK decays 



Measurements of time-dependent CP violation in b — > qqs (q = u, d, s) decays offer a method for determining f3 
alternative to the one based on B —> ccK^*'°. Such decays are dominated by b — >• s loop diagrams, and therefore are 
sensitive to possible new physics contributions appearing in the loops of these diagrams. As a result, the effective /? 
(Peg) measured in such decays could differ from f3 measured in B° — > ccK°. Two final states, the CP-even KgKgKg 
and the CP-odd 4>K° S , are particularly suited for a NP search, as /3 e ff for these decays is expected to be very close, 
within a few %, to /3 measured in B decays to charmonium [14H16j . with very little theoretical uncertainties. The 
related decay mode B ± — > §K ± is also a powerful probe to NP, since it is dominated in the SM by a single b — > s 
penguin amplitude, and its direct CP asymmetry is predicted to be small, Acp ~ (0.0 — 4.7)%. 

BABAR has published in [T7] a detailed study of B — > KgKgKg decays, including a time-integrated, CP-averaged 
Dalitz plot analysis - to extract the amplitude and phases of the various resonant and non-resonant contributions 
to the decay - and a measurement of sin 2/3 e ff from the time-dependent CP asymmetry. In the Dalitz plot analysis, 
events are selected in which three Kg decay to charged pion pairs, have a common vertex and form a B meson 
candidate with tties ~ m B an d A_E 0. Continuum qq background is suppressed by means of a neural network NN 
of four event-shape variables. The B° —> KgKgKg decay amplitude is modeled as the sum of a non-resonant term 
and several intermediate two-body amplitudes; the magnitudes and phases of each term are determined, together 
with the signal and background yields, through a maximum likelihood fit to jties, AE, NN, and to the two Dalitz 
plot variables, s m i n and s max , which are the minimum and the maximum, respectively, among the three Mandelstam 
variables sy, where the indices i, j = 1..3 (i ^ j) denote the three Kg. Evidence is found for contributions to 
the final state by f (980)K s , /o(171Q)iTg, / 2 (2010)i^, XcqK° s and non-resonant K° S K S K S , while the fit likelihood 
does not improve significantly when including in the model additional resonances, like the broad /jc(1500) used to 
explain some older data [18] . The global maximum of the likelihood favors large non-resonant and small resonant 
contributions, with significant destructive interference among them. The signal yield obtained from the fit is 200 ±15, 
while the efficiency, obtained from a simulation of signal events where the magnitudes and phases of the intermediate 
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resonances are fixed to the central values returned by the fit, is e — 6.6%, yielding an inclusive branching fraction 
BF(B° -> K° S K° S K Q S ) = (6.19 ± 0.48 ± 0.15 ± 0.12) x 10~ 6 . 

In the CP violation analysis, events are selected in which one neutral J5 tag decays to a (partially reconstructed) 
flavor-specific final state, thus identifying the flavor of the other B (P s ig) as B° or B at the same time, while B s [ g 
decays to either three K% -> tt+tt" (e = 6.6%) or two K° s -> tt+tt" and one K° s 7r°7r° (s = 3.1%). The CP 
asymmetry between the B° and B decays to three Kg has a dependence on the proper time difference At between 
the two B decays of the form [S sin( AmdAt) — C cos(Am c ;At)], where Anid is the mass difference between the two 
neutral B physical eigenstates. In the SM S — — sin 2/3, while C w 0, since the decay amplitude is dominated by a 
single weak phase term and direct CP violation is thus negligible. S and C are extracted from a maximum likelihood 
fit to At, iues, AE and NN for the selected candidates (see Figure [2]). The signal yield is 263 ±21, and the measured 
values of S and C are: 



S = -0.94+°;^( s tat) ±0.06(syst), 
C = -0.17±0.18(stat)±0.24(syst), 



(1) 
(2) 



where the systematic errors are dominated by the At resolution and the fit bias. They are compatible within two 
standard deviations with those measured in tree-dominated modes such as B° — > J/ipK^, as expected in the SM. 
CP conservation (S = C = 0) is excluded at 3.8<r, including systematic uncertainties, thus providing for the first time 
evidence of CP violation in B° — > K°K°K° decays. 
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FIG. 2: Left: Background-subtracted data distribution (points with error bars) and fit (histogram) of the time-dependent CP 
asymmetry for B — > K l ^K^Kg candidates reconstructed by BABAR [17]. Right: two-dimensional scan of —2AlnC, including 
systematic uncertainties, as a function of S and C. The result of the BABAR analysis of B° — > ccK^*' decays is indicated as 
a white ellipse and the physical boundary (S 2 + C 2 < 1) is marked as a gray line. 



BABAR has also published in (TU] a detailed study of the other B decays to three kaons. The decay B° — > K + K~Kg 
is used, by performing a time-dependent analysis, to measure /3 e ff in B° decays to the CP-odd final state 4>Kg and 
also to the (mostly CP-even) other resonant and non- resonant K+ K~ K° s final states. The decay B° K+R-R+ 
is used to determine the direct CP asymmetry in B + — > <j)K + . In order to disentangle the B — > cf>K contribution 
from other resonant or non-resonant contributions, a Dalitz plot analysis is performed for each of the final states 
considered. In addition, to better understand the Dalitz plot structure of these decays and the possible intermediate 
KK resonances, also B + — > KgKgK + decays are studied, since they have a simplified spin structure due to the 
fact that the two Kg mesons in the final state are forbidden (by Bose-Einstein statistics) to be in an odd angular 
momentum configuration. Events with three charged or neutral kaons kinematically consistent with a B decay 
(common production vertex, m^s ~hib, AE w 0) are selected. Continuum background events are suppressed by 
means of a neural network NN of typically five event-shape variables. The signal B decay amplitude is modeled 
as the sum of a non-resonant and of several resonant contributions, whose magnitudes and phases are determined 
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through a maximum likelihood fit to the Dalitz plot variables, together with tties, AE and NN. Only resonances 
that significantly improve the quality of the fit are kept in the nominal signal decay amplitude: in all the three 
channels no evidence is found for the broad /x(1500). The selection efficiency is around 30% for each of the three 
decay modes. The number of selected signal decays, the measured inclusive charmless branching fractions (excluding 
the XcqK contributions) and the CP-violating parameters are summarized in Table |n] Some of the fit results are 
also illustrated in Figures [3] and [|J The CP asymmetry in B + — > <pK + differs from by 2.8a, and is in slight tension 
with SM predictions (0 - 4.7%). The amount of CP violation in P>° -> (j>K% and in other P>° -> K + K~ K% decays is 
in excellent agreement with SM expectations: /3 ff is very close to j3 measured in B decays to charmonium, and C is 
consistent with zero. 

TABLE II: Results of the B -¥ KKK measurements by the BABAR Collaboration [19]. For neutral B decays, the CP 
asymmetry Acp is defined as Acp = —C. 



Decay 




BF (1(T 6 ) 


A C p 


/3 cff 


B+ -> K+K-R+ 


5269 ± 84 


33.4 ±0.5 ±0.9 


(— 1.7±i;| ± 1.4)% (inclusive) 
(12.8 ± 4.4 ± 1.3)% (<t>K + ) 




B° -> K+K~K% 


1579 ± 46 


25.4 ±0.9 ±0.8 


(-5 ±18 ±5)% 
(-28 ±24 ±9)° (/ (980)tfg) 
(-2±9±3)° (Other) 


(21 ± 6 ± 2)° (4>K%) 
(18±6±4)° (/ (980)Jfg) 
(20.3 ± 4.3 ± 1.2)° (Other) 


B+ -> K° S K%K + 


632 ± 28 


10.1 ±0.5 ±0.3 


(4±5±2)% 
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FIG. 3: Distributions of m ES for B+ K + K~K+ (left), B° -> K+K~K%, and B+ ->■ K%K° S K+ (right) selected by 
BABAR [19] (points with error bars). The fit results are overlaid. 
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FIG. 4: Left: Background-subtracted distributions of B — > K + K~ K ± candidates in data, plotted separately for B + and 
B~ events, in the "(/)(1020) region" where m K + K ~ w (considering the kaon pair with lowest invariant mass). Center and 
right: time-dependent CP asymmetry, as a function of At, for B° — > K + K~ Kg candidates in the 0(1020) region (center) and 
in the </>(1020)-excluded region (right). The points represent signal-weighted data, and the line is the fit model. 
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C. Search for CP violation in three-body B — > K decays 



Recent measurements of rates and asymmetries in B — > Kit decays may hint to possible new physics contribu- 
tions, but hadronic uncertainties prevent a clear conclusion in this direction. Additional information can be obtained 
through a data-driven approach involving measurements of all the observables in the B — ¥ Kir system, or looking 
at the related pseudoscalar-vector decays B — > K*tt and B — » Kp. In such decays, the ratios of tree-to-penguin 
amplitudes are predicted to be 2 to 3 times larger than those in B — > Kir, possibly leading to larger CP asymme- 
tries [2QH22]- 

BABAR has published in [53] a measurement of the branching fraction and CP asymmetry of the decay B ± — > 
K*{892) ± tt , based on the final T(4S) sample. The K*{892)+ is reconstructed in the K+tt° decay channel, and 
discrimination between B + — > K* (892) + tt° decays and other (resonant or non-resonant) B + —> K + tt it decays is 
achieved through inspection of the two-body invariant mass distributions of the particles in the final state, since 
the Dalitz plot of this decay is characterized by rather narrow resonances (K* (892) + , / (980), Xco)- A total B + —> 




FIG. 5: Left: jties distribution of B — > K candidates selected in data by BABAR [23] . after applying a requirement 

on the neural-network variable that retains 60% of the signal while suppressing a large fraction of the qq background. The 
solid (blue) lines represents the tubs projection of the total fit, the dashed (green) lines the total background contribution, and 
the dotted (red) lines the qq component. The dash-dotted lines represent the signal contribution. Right: efficiency-corrected, 
background-subtracted distribution of the invariant mass m K+n a for B — ¥ K tt°tt candidates where m K+7T o is close to 
the K*(892) + mass. m I( + n o . is the mass of the K + ty° combination with lower invariant mass. The vertical red arrows define 
the limits of the signal region while the blue horizontal arrows define the sideband background control regions. The curves show 
the results of the fit used to cross-check the procedure, for the total (blue solid line) and background-only (red dashed line) 
components. 



K + TT°Tr yield of 1220 ±85 events (significance > IOct) is obtained from a ML fit to toes and to a neural network NN 
of event-shape variables (see Figure [5) left). The contribution from B+ K* (892)+ 7T° is obtained by subtracting 
non-_B + — > K + ir°ir° events using the sPlot technique and then studying the distribution of m K +^o , the invariant 
mass of the K + tt° combination with lower invariant mass. The sidebands of mx+^o are used to estimate the non- 
-ftT*(892) + 7r° background in the signal region near the K*{892) + mass, which is then subtracted from the number of 
candidates observed in the signal m K +^o region (see Figure |5j right). The final results are summarized in Table III 
For B ± — > K* (892) ± 7T°, a branching fraction of (8.2 ± 1.5 ± 1.1) x 10~ 6 is measured, with a CP asymmetry close to 
zero (—6%) but also consistent, within errors, with large CP violation (Aqp ~ 20 — 30%) as expected in the SM. 
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TABLE III: Summary of measurements of branching fractions (averaged over charge conjugate states) and CP asymmetries 
measured by BABAR in B + —¥ K + tt°ty° decays [25]. Both product branching fractions and those corrected for secondary 
decays are shown. For each result, the first uncertainty is statistical, the second is systematic and the third, where quoted, is 
the error on \c0 tv°tt°. The notation Rh refers, where applicable, to the intermediate state of a resonance and a bachelor 
hadron. 



Decay 


BF(B+ -> Rh^r K + n°n°)(10- 6 ) 


BF(B + -¥ Rh)(lQ- 6 ) 


A C p 


B+ -> AT+ttV 


16.2 ± 1.2 ± 1.5 




-0.06 ±0.06 ±0.04 


B + -> K*(892)ir° 
B+ -> f (98Q)K+ 
B+ -> Xc oK+ 


2.7 ±0.5 ±0.4 

2.8 ±0.6 ±0.5 
0.51 ±0.22 ±0.09 


8.2 ± 1.5 ± 1.1 
180 ± 80 ± 30 ± 10 


-0.06 ±0.24 ±0.04 
0.18 ±0.18 ±0.04 
-0.96 ±0.37 ±0.04 



IV. SEARCH FOR NEW PHYSICS IN VERY RARE B DECAYS 

Some B meson decays are expected in the Standard Model to have branching fractions below the experimental 
sensitivity of the £?-factories. The observation of such very rare decays would therefore provide evidence for new 
physics beyond the SM. This is the case for instance of lepton-fiavor violating B decays and leptonic B — > vv decays. 

A. Search for lepton flavor violation in B -> hrt 

Lepton flavor violating decays of B mesons can occur in the SM through loop processes that involve neutrino 
mixing. These are highly suppressed by powers of m 2 , /m 2 ^ , and have predicted branching fractions many orders of 
magnitude below the current experimental sensitivity. However, in many extensions of the SM, B decays involving 
lepton flavor violation are greatly enhanced. In some cases, decays involving the second and third generations of 
quarks and leptons are particularly sensitive to physics beyond the SM [24H27] , 

BABAR has presented at the FPCP 2012 conference, and later published in [55], a search for the eight lepton 
flavor violating decays B ± — > h ± r£ , h = tv, K, I = e, fj,, using the final T(45) dataset. Events in which one B 
meson (B tag ) is fully reconstructed in one of several hadronic final states, B tag — > £)W° n 2 K ± n 3 Kg 71471"°, 

are selected. The number of tracks from the signal B (B s i g ) decay and their particle identification information are 
required to be consistent with that of a B ± — > h ± r£ decay followed by a one-prong t decay (producing either a 
charged electron, muon, or pion), in order to improve the signal-to-background ratio. Continuum background is 
suppressed using a likelihood ratio based on event shape information, unassociated calorimeter clusters, and the 
quality of muon identification for channels that have a muon in the final state. Background from semileptonic B or 
D decays is suppressed by rejecting events where two oppositely-charged B s i g daughters are found to be kinematically 
compatible with originating from a charm decay, by computing their invariant mass in the hypothesis of a Kir pair. 
Using the momenta of the reconstructed B, h, and £ candidates, the four momentum of the r lepton in the center- 
of-mass frame is fully determined, thus circumventing the problem of measuring the undetected neutrino(s) from the 
t decay: 

P T = -Ptag -Ph-Pe (3) 

77* JTi* 77'* 171* / A\ 

The resulting r candidate mass, m T = \l E* 2 — \p* T \ 2 , is then used to discriminate against combinatorial background 
(see Figure|6|. A broad m T sideband ([0, 3.5] GeV/c 2 ) is used to estimate the background in the signal region of m T 
close to the PDG r mass, using a simulation of background events to relate the number of background candidates in 
the sideband region to the number of background candidates in the signal region. To reduce systematic uncertainties, 
the signal branching fraction is determined by using the ratio of the number of B — > hr£ signal candidates to the 
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FIG. 6: Observed distributions of the r invariant mass for the B + — > K + t~ u + candidates (left) and the B + — > K + r + e~ 
candidates (right) selected by BABAR [28]. The distributions show the sum of the three r channels (e, fi, n). The points with 
error bars are the data. The solid line is the background MC which has been normalized to the area of the data distribution. 
The dashed vertical lines indicate the m T signal window range. The inset shows the m T distribution for signal MC. 



yield of control samples of B — > D^*'°i v, D° — s- K n + decays in events with a fully reconstructed hadronic B 



tag- 

tag 
Dfi 
tag 



AT & 

BF hTl - BF Dlv — (5) 
^Div e hri e hTl 



Here e^S aim £ £>fi/ are the (similar) B t3bg reconstruction efficiencies for events containing the different B s i g decays, 
while effL and s s ^ gl> are the signal selection efficiencies. A charged Stag candidate is properly reconstructed in 
approximately 0.25% of all BB events. The B — > hr£ selection efficiency varies between 2% and 10% depending on 
the final state; the D^°£v signal efficiencies are close to 50%. The D°£v and D*°iv yields are determined through 
a fit to the energy difference of the selected candidates, 

kE DtlJ = E* K + E* + Eg + El - £b cam , ( 6 ) 

where E* = p* y = \ — p£ ag — p* K — — | . The energy difference is peaked at for B — > D°li> and is peaked at lower 
values (« -200 MeV) for B D* Q tv (because of the undetected 7 or tt° from the D*° decay), while B D**°£i> 
and other backgrounds do not peak in AEuiv (see Figure [7]) . 

No evidence for the B — > hr£ decays is found, and 90% confidence level upper limits are set - with a frequentist 
technique - on each branching fraction at the level of a few times 10~ 5 , as summarized in Tabic fV The results are 



used to improve model-independent limits on the energy scale of new physics in flavor-changing sermonic effective 
operators to A bd > If TeV and A bs >15 TeV. 



B. Search for new physics in B decays to invisible final states and to uwy 

The decay B a — > vv, which would give an invisible experimental signature, is strongly helicity-suppressed in 
the SM by a factor of order (m^/ms) 2 and its branching fraction is thus well below the B factory experimental 
observability [29] ■ The SM expectation for the B° — > i/Dj branching fraction is predicted to be of order I0~ 9 , with 
very little uncertainty from hadronic interactions [30 . An experimental observation of an invisible (+7) decay of a 
B° at the B factories would thus be a clear sign of physics beyond the SM. Supersymmetric models or models with 
large extra dimensions allow significant, although small, rates for invisible B° decays, with branching fractions up to 
the I0- 7 - I0" 6 range [HHSSl- 
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FIG. 7: Distributions of the three-component AEoe^ unbinned maximum likelihood fits of the data for the B — > D^*'°ev (left) 
and B — > D^*' ^iv (right) control samples selected in the B —¥ hr£ BABAR analysis [2S]. In each plot, the points represent 
the data, the solid blue curve is the fit result, the long-dashed green curve is the D° component, the dash-dotted purple 
curve is the D*° component, and the dotted blue curve is the component, which also includes any residual combinatorial 
background. 



TABLE IV: Branching fraction central values and 90% C.L. upper limits (UL) for BF(B+ — > h+r€). The uncertainties 
include statistical and systematic sources. 



Mode 


BF(B -» h- 
central value 


(xHT 5 ) 
90% C.L. UL 


B+ -> K + t~h+ 


o.8 tli 


< 4.5 


B + -> K+r+pT 


-o.4 ±S:I 


< 2.8 


B + -> K+r~e+ 


0.2 t?;J 


< 4.3 


B + -> K+T+e' 


-i-3 


< 1.5 


B+ -+-K+T-U+ 


0-4 ±1.1 


< 6.2 


B + -> 7T+T+M" 


o.o ±2:g 


< 4.5 


B+ -> 7T + T~e+ 


2-8 tn 


< 7.4 


B+ -» -K+r+e- 


-3.1 t 2 2 ;l 


< 2.0 



Both Belle and BABAR have presented at the FPCP 2012 conference preliminary results on invisible P>° decays; 
BABAR has also searched for B° —> decays. The results have been later publicly documented in [M] and [55] . 
The analysis technique is similar for the two experiments: only events where a neutral B tag meson is found are 
considered, and the rest of the event (the "signal side") is checked for consistency with an invisible decay or a decay 
to a single photon of the other neutral B. Events with extra tracks or neutral pions or candidates are vetoed. 
The total energy -E GX tra (or -Eecl) m the electromagnetic calorimeter, computed in the CM frame and not associated 
with neutral particles or charged tracks used in the B tag reconstruction, is used to discriminate between signal 
and background. For B° — > invisible + 7, the energy of the highest-energy photon remaining in the event is also 
removed from the -E e xtra computation. In signal events -E C xtra is strongly peaked at zero, while for the background 
its distribution increases uniformly with -E^tra- The main difference between the Belle and BABAR analyses is that 
Belle performs a full reconstruction of the B tag candidate in hadronic final states, while BABAR opts for a partial 
reconstruction in the D^*'{+v decays. The latter has the advantage of a higher reconstruction efficiency with respect 
to the former (by a factor w 10: the B — > invisible efficiency is 0.18% for BABAR and 0.022% for Belle), but has the 
disadvantage of the presence of the invisible neutrino, which prevents the exploitation of kinematic variables such 
as the reconstructed B° mass. S/B is therefore worse in the case of the partial B tag reconstruction, though the 
background contamination is mitigated by the presence of a high momentum lepton. 
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The signal yield is extracted from a maximum likelihood fit to the two-dimensional distributions of -Eecl and 
of an event-shape variable, cosOb (Belle), or to the one-dimensional -E^tra distribution (BABAR), as illustrated in 
Figures [8] and [9] No evidence for B — > invisible(+j) decays is found, and 90% confidence level upper limits are 
set - with a Bayesian technique - on each branching fraction, at the level of a few times 10 -5 , as summarized in 
Table [V] The Belle limit is a factor m 5 worse than BABAR because of the lower efficiency of the hadronic tag full 
reconstruction. 




FIG. 8: -Eecl distribution of B — > invisible candidates (points) selected by Belle, with fit results superimposed. The red 
cross-hatched region is the signal component, on the top of the total background shown in the yellow filled histogram. The 
blue dashed curve is the generic B contribution, which is larger than the total because of the negative fit result for the non-B 
background shown in the green dotted histogram. The purple hatched area corresponds to the rare B contribution. 




FIG. 9: -Ecxtra distributions of B — > invisible (left) and B — > invisible + 7 (right) candidates selected by BABAR (points), 
with fit results superimposed: signal (purple dotted line), background (green, dotted line) and sum of signal and background 
(blue solid line). 



TABbE V: 90% C.L. upper limits on the branching fractions of B — > invisible(+ r y) decays measured by the Belle and BABAR 
collaborations [34l 135] . 



Decay 


BF (Belle) 


BF (BaBar) 


B° — > invisible 


< 13 x 10 -5 


< 2.4 x 1(T 5 


B° — > invisible + 7 




< 1.7 x 1(T 5 
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V. CONCLUSION 

Several results on rare B decays obtained at the B factories have been discussed. The main results are: 

• CP violation has been observed with > 3<r significance in B —¥ r\K ± and B ± — ¥ rpr . Large CP asymmetries 
are observed. B°^r)K° is observed for the first time, with a branching fraction smaller than that of B + ^nK + . 

• CP-violation parameters in B a — ¥ KgKgKg are compatible within two standard deviations with those measured 
in tree-dominated modes such as B° — ¥ J/tpKg, as expected in the SM. CP conservation is excluded for the 
first time in these decay with a significance of more than 3cr, including systematic uncertainties. 

• The CP asymmetry in B + — ¥ cj>K + differs from by 2.8a, and is in slight tension with SM predictions (0—4.7%). 

• The amount of CP violation in B° — ¥ 4>K° S and in other B° — > K + K~Kg decays is in excellent agreement with 
SM expectations: /3 e ff is very close to /3 measured in B decays to charmonium, and C is consistent with zero. 

• No evidence of CP violation has been found in B ± —¥ X* (892) ± 7r°, though the large uncertainty on the CP 
asymmetry does not exclude CP-violating effects as large as 20 — 30%, as expected in the SM. 

• No evidence of lepton-flavor violating B ± — ¥ hr^rt decays has been found, and upper limits of a few times 
10 -5 on their branching fractions have been set, improving by a factor 5 previous limits on \x — r lepton-flavor 
violating effective couplings. 

• No evidence of B — ¥ invisible(+j) decays has been found, and upper limits of a few times 10~ J on their 
branching fractions have been set, 1-2 order of magnitude still higher than the values predicted by some 
theoretical models beyond the SM. 

All the results are in agreement with Standard Model expectations: no evidence of new physics is found. 
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